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Double open-charm production is one of the most promising channels to disentangle single from 
double parton scattering (DPS) and study different properties of DPS. Several studies of the DPS 
contributions have been made. A missing ingredient so far has been the study of polarization 
effects, arising from spin correlations between the two partons inside an unpolarized proton. We 
investigate the impact polarization has on the double open-charm cross section. We show that the 
longitudinally polarized gluons can give significant contributions to the cross section, but for most of 
the considered kinematic region only have a moderate effect on the shape. We compare our findings 
to the LHCb data in the D°D° final state, identify observables where polarization does have an 
impact on the distribution of the final state particles, and suggest measurements which could lead 
to first experimental indications of, or limits on, polarization in DPS. 


I. INTRODUCTION 

Processes in hadron collisions where two partons from each hadron take part in separate partonic subprocesses, 
double parton scattering (DPS), contribute to several final states of interest at the LHC. DPS is a relevant background 
to precise Higgs boson coupling measurements and searches for new physics [1-4]. The theory for DPS is still frag¬ 
mentary, but major improvements have been made over the last couple of years moving towards a reliable description 
within perturbative QCD [5-8]. Despite this development there are still several important questions which have to 
be worked out. 

For sufficiently inclusive cross sections DPS is formally a power suppressed contribution, but in certain regions of 
phase space double and single parton scattering contribute at the same power [9]. Even for inclusive cross sections, 
DPS can in specific situations compete with single parton scattering - for example when the single parton scattering 
is suppressed by multiple small coupling constants. DPS is increasingly relevant with collider energy, and will hence 
be further enhanced when the LHC restarts to collide protons at larger center of mass energies. The reason is the 
rapid increase of the density of partons with energy and towards smaller x-fractions. 

DPS signals have been measured at the LHC by both ATLAS [10] and CMS [11] in the IV-boson plus dijet final 
state. Of particular interest for our present study is the LHCb measurements of double open-charm production [12], 
in final states such as D°D°. Among the most promising channels for a clean separation of double from single parton 
scattering are the production of two same sign Pb-bosons and double open charm quarks [13-19]. In fact, studies have 
shown that for double open-charm production in the kinematical region of the LHCb measurement, double parton 
scattering dominates over single parton scattering [17]. 

DPS cross sections are factorized into two hard partonic subprocesses and two double parton distributions (DPDs). 
Only little is known about the size of the DPDs. They have been studied in a variety of quark models [20-23], 
including correlations between the two partons inside the same proton. The correlations have generally been found 
to be sizable. An open question in double parton scattering is the effects of quantum-number correlations. These 
include correlations between the spins, colors, flavors and fermion numbers of the partons [5, 9, 24, 25]. Upper limits 
on the DPDs describing quantum-number correlations have been derived [24, 26]. For polarized DPDs these limits 
have further been shown to hold under radiative corrections from the leading-order double DGLAP evolution up to 
higher scales. 

In particular the spin correlations (described by polarized DPDs) have direct relations to the directions of the final 
state particles, and thus have the potential to change both the sizes of the DPS cross sections and the distributions of 
the produced particles. For example, azimuthal modulations have been found for double vector boson production [27]. 
The effects of the quantum correlations on DPS cross sections have been calculated [5, 27] but so far no numerical 
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results at the cross section level have been obtained. Through studies of the scale evolution of the DPDs, limits on 
the degree of polarization and thereby its possible effect on DPS cross sections at different scales were set in [26]. 

In this paper we examine the effect that polarization in DPS can have on the double cc production in kinematic 
regions resembling those of the LHCb D°D° measurement [12]. Several studies of this process already exist in the 
literature, but so far all have neglected the possibility of spin correlations. We demonstrate for the first time the 
quantitative impact of polarization on any DPS cross section. 

The structure of the paper is as follows: In section II we discuss some basics of DPS with focus on polarization, 
introduce the different polarized and unpolarized double gluon distributions and discuss their scale evolution. In 
section III we present the analytical results for the cross section calculation including all possible polarizations of two 
gluons in an unpolarized proton. In section IV we discuss the models for the DPDs which we use in order to obtain 
numerical results - which we present and compare to LHCb data in section V. We summarize our findings and discuss 
their implications in section VI. 


II. DOUBLE GLUON DISTRIBUTIONS 

Under the assumption of factorization, as illustrated in figure 1, the DPS cross section can be expressed schematically 
as 


da 


IL=i dxidxi 


DPS 


= ^2 


d 2 y F(x 1 ,x 2 ,y)F(x 1 ,x 2 ,y ), 


(1) 


where di represents hard subprocess i and C is a combinatorial factor equal to two (one) if the final states of the two 
subprocesses are (not) identical. F (F) labels the double parton distribution of the proton with momentum p (p). 
The DPDs depend on the longitudinal momentum fractions of the two partons Xi ( Xi ) and their transverse separation 
y. No complete proof for factorization in DPS exists, but several important ingredients have been established [5, 9]. 
The cross section expression (1) is schematic as the labels for the different flavors, colors, fermion numbers and spins 
of the four partons are implicit. The possibility of interference between the two hard interactions, and correlations 
between the two partons inside each proton renders this structure significantly more complicated in DPS than for the 
case with only one hard interaction. Of particular interest for our purposes are the correlations between the spins of 
two gluons, and between the spins and the transverse separation, which lead to polarized gluon DPDs. 



Figure 1. Gluon induced double parton scattering. The green fields represents the two DPDs while the blue and red fields 
represent the first and second hard interaction. Xi ( Xi ) are the longitudinal momentum fractions of the partons from the proton 
with momentum p (p). y and 0 are the transverse positions of the first and second hard interaction. 

For the DPDs describing two gluons in an unpolarized right-moving proton we write [9] 

F aia2 ( Xl ,x 2 ,y) = 2p+(x 1 p + )- 1 (X 2 P+)- 1 J dy~ e ^+^)p + 

x {p\Oa 2 {0,z 2 )O ai (y,zi)\p) . (2) 

We use light-cone coordinates v± = (v° ± v 3 )/\/2, bold font to denote the transverse component v = (v 1 ^ 2 ) of any 
four-vector v and vt = |v|. The operators expressed in terms of the gluon field strength tensor reads 

Oat (y, Zi) = njf G + i' (y - \ Zl ) G+i (y + §*) 


zf =y + = 0, z i =0 


( 3 ) 
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with projections 


nf =8”\ = ie jjl , [^Sg] jj ' = T jj,,kk ' (4) 

onto unpolarized gluons (g), longitudinally polarized gluons (A g) and linearly polarized gluons (Sg). The tensor 

T jj '’ kk ' = i (S jk S j ' k ' + 6 jk '6 j ' k - S jj '6 kk ') (5) 

satisfies t- 7 - 7 ,kk r kk ’ 11 = t- 7 - 7 ’ 11 and is symmetric and traceless in each of the index pairs (jj r ) and (kk'). 

A decomposition of the nonzero distributions for two gluons in terms of real-valued scalar functions has been given 
in [26] 


F gg (x 1 ,x 2 ,y) = f gg {x 1 ,X 2 ,y), 
FAgAg(Xl,X 2 ,y) = fAgAgixi, X 2 ,y) , 

F g6g{ x l’ x 2>y) = T n '’ VV M‘ 2 fgSg{x 1 ,X2,y) , 

F Sgg(x 1 ,X2,y) = T n '' VV M 2 fsgg(Xl,X 2 i y) , 

F SgSg k ( X 1,X 2 ,y) = ^T n ’’ kk ' f S gSg(xi,X 2 ,y) + 


( T n',yy T kk',vy _ T n',yy T *k',yy) M 4 f t SgSg (x 1 ,x 2 ,y ), 


( 6 ) 


where M is the proton mass and y J = e n y J . The notation where vectors y or y appear as an index of r denote 
contraction, i.e. t- 7 - 7 vv = r- 7 - 7 ,fcfc y k y k etc. The distributions of longitudinally polarized gluons carry open transverse 
indices j,j\k,k' = {1,2} corresponding to the polarization vectors of the gluons which are contracted with the 
partonic cross sections. 

The double parton distribution f gg represents the probability of finding two gluons with momentum fractions X\ 
and x 2 at a transverse separation y. The distribution of longitudinally polarized gluons fAgAg describe the difference 
in probability between finding the two gluons with their helicities aligned rather than anti-aligned, while linearly 
polarized gluons are described by hclicity interference distributions, see for example [26, 28] in the context of DPS. 


A. Evolution of the double gluon distributions 


The scale evolution of the DPDs is governed by a generalization of the DGLAP evolution equations. Two versions 
exist in the literature: one homogenous equation describing two independent branchings of the two partons, and 
another including the splitting of a parent parton into the two partons which subsequently undergo hard scatterings 
[29-33]. Which one is the correct one for describing DPS is still under debate [6, 9, 34-42]. The contribution from 
the splitting term was investigated in [15] for double cc production and was seen to give a sizable contribution to the 
cross section, but also that the perturbative splitting preferred to take place early on - and evolve as two separate 
branches for most of the evolution range. Including such a term in our study could naturally lead to an enhancement 
of the effect of the polarization and we will return to this discussion in section VI. In the following we will make use 
of the homogeneous version, under the assumption that the physics of the single parton splitting contribution can be 
treated separately. 

The evolution equation for the unpolarized double gluon distribution then reads 


df gg {xi,x 2l y, y) 
d In /i 2 


2 n ^ ' [ F ga ®1 fag + Pga ®2 fga] > 

a=g,q,q 


(7) 


where 


P a b{ ■) ®i /be ( -,x 2 ,y;y) 
Pab{ ■) ®2 fcb{xi, - ,y\y) 



X2 du l n 

( X\ 

- Pab 

— 

Ul 1 

{Ul 

K1 du 2 „ 

( X2 

- Pab 

— 

u 2 

{U2 


fbc{ui,x 2 ,y;y ), 

fbc(xi,u 2 ,y,y) 


( 8 ) 


are convolutions in the first and second argument of the DPDs with the leading-order splitting kernels P a b known 
from DGLAP evolution of single-parton distributions. Polarized DPDs follow equivalent evolution equations with the 
splitting kernels replaced by their polarized analogues. A more thorough discussion of the evolution of the polarized 
DPDs and expressions for all splitting kernels are given in [26]. The evolution of the unpolarized gluon distribution 



4 


leads to a violent increase at low momentum fractions, in particular at low scales where the QCD coupling constant 
is large. This is due to the 1/x behavior of the unpolarized splitting kernel in the limit where x tends to zero. The 
splitting kernel for a longitudinally polarized gluon on the other hand approaches a constant in this limit, while the 
one for linearly polarized gluons goes as x. The polarized distributions therefore do not experience this rapid increase 
and evolution will suppress the relevance of polarized gluons - in particular the linearly polarized ones. The rate at 
which this suppression takes effect leads to the expectation that at large scales (and not too large x) polarized gluons 
can be neglected in phenomenological calculations of DPS cross sections [43]. However, for double cc production the 
scales are low and there is only little room for evolution. This motivates the study of the effects of polarization in 
this process, and could, when confronted with experimental results lead to the first measurements of, or limits on, 
polarization effects in DPS. 


III. DOUBLE cc CROSS SECTIONS 


We next present the analytic results of the cross section calculation, dividing the results into contributions from 
the different polarizations. The non-zero results come from gluons which are unpolarized, longitudinally polarized, 
mixed unpolarized-lincarly polarized and purely linearly polarized. 

Following experimental conventions, we present our results in the center-of-mass (CM) frame of the two protons, 
with 2 -axis along the proton with momentum p and i-axis as pointing towards the centre of the LHC ring. With 
this choice of i-axis, without reference to any direction defined by the process itself, any azimuthal dependence must 
show up as differences between the azimuthal angles describing the transverse directions of the final state particles. 
Double cc production in the kinematic region of interest is, to good approximation, initiated by gluons (see e.g. [17]). 
Therefore we limit ourselves to the partonic subprocesses of figure 2, where the cc systems are produced by s-channel 
gluons or t-channel (u-chanel) charm quarks. 



Figure 2. Leading order Feynman diagrams contributing to each of the two partonic subprocesses. 


The cross section contribution from unpolarized gluons reads 
d a (gg)(gg) 


dyidy 2 d 2 p 1 d 2 p s 


1 (o 2 2 N t 
= - 87T a 


N? - 1 


dxidx 2 


1 


2XjXi 


167 r 2 sf 2Xi — XTiG Vi 


(1 - + zl - l/N 2 

(1 - Zl)Z! 

(1 - z 2 ) 2 + 4 - l/N 2 
(1 - 22)^2 


(1_ z 2 ) 2 + 2 2 + 4zi( i_ Z i) ( 1 _^\ 

' TTbrj-i-t / 


TTtrjp-^ / TYl/rj-i-^ 


(1 — 2 2 ) + z 2 + 422(1 — 22 ) I 1 ~ 


9 \ 9 

m \ m 


t 2 / " b T 2 J 


x / d 2 yfgg(x 1 ,x 2 ,y)f gg (x 1 ,x 2 ,y), 


(9) 


where yi and p t are the rapidity and transverse momentum of the charm quark produced in interaction i = 1,2. The 
variables in the cross section are given by 


— 




XjXTie Vi 
2Xi — XTie Vi ’ 


XTi = 


2rriTi 


mn = \ P 2 + m 2 


( 10 ) 


§i and ti are the usual Mandelstam variables of the partonic cross section i. s is the center of mass energy of the 
proton collision, m is the charm mass and y, is the rapidity of the charm quark from interaction i. 
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For gluons with longitudinal polarization the cross section is 


^ cr (AgAg)(AgAg) 

dyidy 2 d 2 p l d 2 p 2 


1 

2 

x 

x 


87 r 2 ai 


N r 


Nq 1 
(1 - Zl ) 2 + z 2 - 

(1 - Zi)zi 
(1 -Z 2 ) 2 + Z 2 - 

(1 - z 2 )z 2 



, 

rr 

1 

2 XiXi 

dx ^ 

' 167r 2 s 2 : 

i=l L 

2 Xi — XTi& Vi 

1 - 

2TO 2 N 
2 

IK 1 - 

-*l) 2 +^] 


m Tl / 


1 - 

2m 2 N 

m T2 J 

) K 1 - 

- Z 2) 2 + zl] 


X 


J d 2 yfAgAg(xl , X 2 , y)fAgAg {x\ ,X 2 i y). 


( 11 ) 


Worth noticing is that differences in the partonic cross section between longitudinally and unpolarized gluons are 
suppressed by m 2 /??!^. For transverse momenta of the outgoing charm quarks above a few GeV, this suppression is 
strong and already at pn = 3 GeV we have m 2 /m‘^ i = 0.16. This tells us that differences in the distributions of the 
final state charm quarks produced at large pr between longitudinally polarized and unpolarized gluons will to good 
approximation originate in differences between f gg and fA g A g - 

Unpolarized mixed with linearly polarized gluons gives the cross section 


d <J (Sgg)(g$g) 

dy 1 dy 2 d 2 p 1 d 2 p 2 


2„2 


8n a' 


N r 


dx\dx 2 


1 


2 XiXi 


m 


s N 2 1 

x ((1 — z\) 2 + zf — 1/-/V 2 ) 
x ((1 ^ z 2 ) 2 + z 2 - 1/N 2 ) 


i=l 

2 


167t 2 s 2 2xi — XTi& Vi 
2 


T1 


1 - 


TO 


"T 1 


TO 


T 2 


1 - 


TO 


T 2 


x cos(2A (j)) / dyy M fs gg {xi,x 2 ,y)f gSg (xi,x 2 ,y), 


( 12 ) 


and the same result for the case when the linear and unpolarized gluons are interchanged, i.e. with g o 5g. A <fi = 
<f>i — </> 2 , where <f>i is the azimuthal angle of the outgoing c-quark from hard interaction i. The cos(2A </)) dependence is 
an effect of the difference in helicity between the amplitude and conjugate amplitude for the linearly polarized gluons. 
This term gives rise to the same kind of modulation in the azimuthal angle as observed by LHCb in the D°D° final 
state [12]. However, in our leading order cross section, the whole contribution for the mixed linear-unpolarized gluons 
is suppressed by m 2 /m^ i for each of the two hard subprocesses. This suppression gives low analyzing power, and 
indicates already at this level that the contribution should be small. The suppression arises in the terms where there 
is an helicity flip in the hard cross section. For zero quark masses these terms in the partonic cross section tend to 
zero. The nonzero charm mass allows for a nonzero result of the mixed un- and linearly-polarized gluons, but only 
at the price of the suppression factor. This has previously been discussed in the context of heavy quark production 
with transverse momentum dependent parton distributions [44] . It is however interesting to note that the suppression 
could be lifted if the gluons are given a transverse momentum, for example by radiating off a gluon. The next-to- 
leading-order (NLO) correction to the cross section is expected to be large [17], and a large NLO contribution in 
combination with a lifting of the suppression have the potential to result in a large enhancement of this contribution. 
We will return to this point in the discussion of section VI. 

The cross section for gluons with linear polarization is 


d a (5g5g)(5g5g) 

dy 1 dy 2 d 2 p 1 d 2 p 2 


— 87t a 


2„2 


IV, 


IV?- 1 


dx\dx 2 


2 XiXi 


i =1 


167 t 2 s ? 2 Xi — XTie Vi 


x ((1 - Zl ) 2 + zl- 1/N 2 ) ((1 - z 2 ) 2 + z 2 - 1/N 2 ) 


1 - 


°T1 


1 - 


m 


m; 


t 2 


cos(4A</)) 


m 


m. 


T1" l T2 


x 


J d 2 yfsg8g(x 1 ,x 2 ,y)f8 g 8g{x 1 ,x 2 , y). 


(13) 


The pure linearly polarized contribution to the cross section has a cos(4A</>) dependent part and a A</> independent 
part. The A<^> independent term comes with a double helicity flip in the two hard cross sections and is heavily 
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suppressed. Notice that the fg gSg term in (6) does not contribute. This is because it results in a dependence on the 
angle between the directions of the outgoing charm quarks and the direction y between the two hard subprocesses, 
which vanishes upon integration over y. 


IV. SIMPLE MODEL FOR DPDS 

In order to obtain numerical results we need an initial ansatz for the DPDs at some low starting scale. We decompose 
the unpolarized DPDs into two single parton distributions and a y dependent function assumed to be universal and 
Xi independent 

fgg(xi,X 2 ,y,Qo) = f g (xi,Qo)f g (x 2 ;Qo)G(y). (14) 

This is an ansatz commonly used for DPS phenomenology but its validity is questionable and in some kinematic 
regions wrong. The easiest way to see this is in the region of large x^. Momentum conservations forces x\ + x 2 < 1 
on the left side of (14), but the right hand side does, as it stands, not respect this constraint and gives nonzero values 
as long as both momentum fractions individually are below 1. A way to reinstate this limit is to multiply the ansatz 
with the factor (1 — X\ — x 2 ) to some positive power, however for the charm production the contribution of the large 
Xi region is negligible and we apply a strict cutoff at the kinematic limit. Despite its limitations, the ansatz provides 
a useful starting point for DPS studies and we use it as input for the unpolarized DPDs at some low starting scale. 
These input distributions will then be evolved to higher scales with the double DGLAP equations. The numerical 
results will only be given in terms of ratios of cross sections, in which the y dependence cancels. For unpolarized 
distributions, the difference between separately evolving the two parton distribution functions (PDFs) or evolving the 
DPD with the factorized initial ansatz is small, except in the large Xi region [43]. 

For polarized distributions, which describe the correlation between the spin of the two partons, it does not make 
sense to decompose it into polarized single parton distributions - which describe the correlation between the spin of 
one parton and the spin of the proton. Instead we use the positivity bounds in [26] to set upper limits on the sizes of 
the polarized distributions in terms of the unpolarized. We are interested in examining the maximal effects possible 
from the different polarizations and therefore saturate the bounds for each polarized DPD independently. This results 
in the relations 

fAgAg{X\,X 2 ,y, Qo) = f gg {x\, x 2 , y , Qo), 
fgSg{xi,x 2 ,y 1 Qo ) = y 2 M 2 f gg (x 1 , x 2 , y; Q 0 ), 

fs g 5 g{xi,x 2l y,Q 0 ) = y 4 M 4 f gg (xx , x 2 , y ; Q 0 ) (15) 

at the initial scale Qo ■ If the bounds are fulfilled at an initial scale they will remain valid at all larger scales, but 
typically be violated at lower scales. We therefore choose to saturate the bounds at some low Qo, and use the double 
DGLAP evolution (7) with polarized splitting kernels to obtain the polarized DPDs at higher scales. 

A larger Qo gives less room for evolution and therefore less suppression of the polarized contribution to the cross 
section. Qo should be chosen such that one is in a regime where perturbative QCD is expected to give sensible results. 
For the usual PDFs, the starting scale is often chosen somewhere around 1-2 GeV, and several of the leading order 
distributions go negative when evolved backwards to scales below 1 GeV. Another issue is the large uncertainty even 
for the distribution of a single unpolarized gluon at small scales and momentum fractions. The smaller values we take 
for the initial scale the larger this uncertainty is, and we would like to stay at a scale where we can compare between 
different sets of single parton distributions. We use this as a guidance in choosing starting scales at which to saturate 
our polarized bounds, and conclude that a choice somewhere between 1 and 2 GeV is reasonable. We investigate the 
impact of this choice by varying the input scale between the two values. 

For input PDFs we will use the leading-order GJR distributions [45]. At these scales there are still large differences 
between different PDF sets, and we have investigated how they influence our results by switching to the MSTW20081o 
distributions [46]. There are clear differences, especially with Qo — 1 GeV, but the differences are smaller than those 
obtained by changing from Qo = 1 GeV to Qo = 2 GeV. As a rule of thumb, the MSTW distributions give smaller 
polarization than the GJR distributions - see [43] for more details on the effects of changing between different sets of 
PDFs. We use the values of a s and the charm mass m used in the PDF sets for concistency. 

In addition, we examine the effect of changing our modeling of the polarized DPDs. Instead of taking the maximal 
allowed polarization, we can create a model built on ratios of splitting kernels describing the branching of a parent 
parton into the two gluons which subsequently undergo hard scatterings. For the longitudinally polarized DPD this 
results in 

_ Tg^AgAg J. _ ZZ (2 — ZZ ) f 

^ Tg^gg ^ _ Z^ + Z^ + Z W " ^ ' 
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where z = X\j{x\ +X2) and z' = 1 — z. This model has been described in more detail in [43], which gives a complete list 
of expressions for the different DPDs. We will only display results obtained in this model for longitudinally polarized 
gluons. 


V. NUMERICAL RESULTS AND COMPARISON WITH DATA 

We next turn to the numerical evaluation of the cross section in the kinematic regions probed by the double 
open-charm measurement by the LHCb Collaboration [12]. The two charm quarks are required to have a transverse 
momentum in the region 3 < pn < 12 GeV and rapidities in the range 2 < y, < 4, at i/s = 7 TeV. The phase-space 
of the two anticharm quarks is integrated over, without experimental cuts, since they remain undetected. 

The DPD evolution equations are evaluated by the code described in [32], which has been modified to suit our 
purposes as described in [43]. The main modifications are the use of the homogeneous evolution equations and the 
incorporation of the polarized splitting kernels for the evolution of the polarized DPDs, listed in appendix A of [26]. 
We generated gridfiles for the DPDs in the range 10 -6 < Xi < 1 with 240 gridpoints in each direction, and 60 points 
in In /i 2 in the range Qq < p 2 < 2 x 10 6 GeV. The phase-space integrations were performed numerically. 

Care must be taken when comparing the data to the results of our calculation. While the calculation produces 
two pairs of charm-anticharm quarks, out of which only the two charm quarks are measured, the data is for D°D°. 
Simply interpreting the variables of the charm quarks as those of the final state mesons neglects the effects from 
hadronization/fragmentation. The assumption that the direction of the charm quark is approximately equal to that 
of the D° is commonly made [18]. For the effect on the absolute size of the transverse momenta, the approximation 
is less accurate, but charm fragmentation functions typically peak around rather large z values [47, 48]. However this 
approximation on the normalized cross section is not likely to change the spectrum at the level of precision we are 
interested in here. Normalizing the results to the total cross section cancels the effects on the absolute size, such as 
the branching ratio of c —> D°. Our primary purpose is not to make exact predictions for the D°D° cross section, 
but rather to examine the effects that polarization has on double charm production. 




3 4 5 6 7 8 9 10 11 12 3 4 5 6 7 8 9 10 11 12 

p T [GeV] p T [GeV] 

Figure 3. Normalized cross section vs the transverse momentum of one of the charm quarks at p = 2m (left) and p = wt 
( right). Overlaid are the LHCb D°D° data [12]. The lower panels show the relative size of the polarized contribution. 


Figure 3 shows the dependence of the normalized cross section on the transverse momentum of one of the two charm 
quarks, pt, as well as the ratio of the polarized over unpolarized contribution. The left panel shows the results with 
the scale choice of p = 2 m and the right with p = mr, where 

m,Ti + rriT2 


rriT = 


2 


(17) 


















is the average transverse mass of the two charm quarks. For both cases, we make visible the dependence of the result 
on the choice of input scale by displaying the results for Q 0 = 1 GeV and Qq = 2 GeV, as discussed in section IV. The 
two lower panels show the relative size of the polarized contribution compared to the unpolarized. The cross section 
result in figure 3 reproduces the data reasonably well. The shape of the cross section only has a tiny dependence on 
the choice of Q o, whilst the contribution from the polarized distributions changes with Qq. Likewise, there is little 
difference in the shape of the cross section with the two scale choices, but the polarized contribution is larger for 
p = 2m than for p = my- This is expected since the latter choice allows for a larger evolution range and thus a 
stronger enhancement of the unpolarized over the polarized DPDs. With p = mr the suppression due to evolution also 
increases with the pr, since increasing pt increases Tot, counteracting the enhancement of the polarized contribution 
from the partonic cross sections. The relative size of the polarized contribution does have a small dependence on pt 
with /i = 2m, but is rather flat for p = mr- 
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Figure 4. Same as in figure 3, but for the cross section differential in the rapidity difference Ay between the two charm quarks. 


Figure 4 shows the dependence of the normalized cross section on the rapidity difference between the charm quarks. 
The cross section results are stable under variations of the scales and nicely reproduce the shape of the data. The 
two input scales have a strong impact on the size of the polarized contributions. The relative polarized contribution 
displays no dependence on the rapidity difference. With p = 2m the ratio of polarized over unpolarized is 30% for 
Qq = 2 GeV and around 4% for Q 0 = 1 GeV. Changing to p = my decreases the ratio to about half, 15% for 
Qo = 2 GeV and 2% for Qq = 1 GeV. 

The cross section dependence on the invariant mass of the two charm quarks, M cc , is shown in figure 5. As in 
the previous figures, the data is rather well reproduced by the double charm cross section calculation. The polarized 
contribution has some dependence on M cc with p = 2m and thus a small impact on the shape of the cross section, 
but this effect disappears for p = mr- With the lower input scale the polarized contribution is a few percent. With 
the larger Q 0 the polarized gluon contribution is 30% of the unpolarized at small M cc and increases up to 50% at 
large M cc for p = 2m, while p = mr gives a ratio just above 10% in the entire M cc range. 

An intriguing aspect of the D°D° results is the azimuthal correlation between the two mesons. This correlation 
differs from that observed between meson final states with an equal number of charm quarks and anti-quarks, such as 
D°D° and D ± D ^ [12], which are dominated by single parton scattering. The angular modulation in D°D a resembles 
that of a cos 2 A (f> dependence, which is naturally produced by the DPS cross section involving a mixture of unpolarized 
and linearly polarized gluons. In our LO calculation, the polarized contribution can still be sizable as demonstrated 
in figure 6, but the large contribution originates from the longitudinally polarized gluons and is thus independent of 
A cj>. 

Instead of the maximal polarization model for the DPDs, we can use the splitting model - where the ratios of the 
perturbative splittings of an unpolarized parent parton into two (unpolarized or polarized) partons are the basis of 
the relations between the different DPDs as explained in section IV. This scenario has a smaller polarization, but 
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Figure 5. Same as in figure 3, but for the cross section differential in the invariant mass of the two charm quarks M cc . 
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Figure 6. Same as in figure 3, but for the cross section differential in the azimuthal angle between the two charm quarks. 


different Xj dependences of the unpolarized compared to polarized DPDs - which could show up as shape differences 
in the rapidity spectrum. Figure 7 shows the cross section as a function of the rapidity difference, with the polarized 
DPDs from the splitting model for Q o = 2 GeV and p, = 2m. The size of the polarized cross section is reduced in the 
splitting scenario to about 5-10%, with some dependence on the rapidity difference. 

Extending the kinematic region to examine in particular the effects of going down towards lower values of pt, 
figure 8 shows the cross section dependence on p t in the kinematic range 1 < pti < 12 GeV. This decreases the 
size of the polarized contribution, which is not surprising since the longitudinally polarized cross section in (11) has 
a (1 — 2 m/mTi) factor for each of the two partonic processes, which decreases when going to smaller pn■ Although 
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Figure 7. Normalized cross section in the splitting scenario vs the difference of rapidity between the two charm quarks at 
p = 2 m (left). Overlaid are the LHCb D°D° data [12]. The lower panel shows the relative size of the polarized contribution. 
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Figure 8. Normalized cross section vs the transverse momentum of one of the charm quarks in the extended region down to 
p Ti = 1 GeV, at /x = 2 m (left) and p = mr (right). The lower panels show the relative size of the polarized contribution. 


the polarization in this region is rather small, the effect on the shape of the cross section is interesting. In figure 9 
we show the double differential cross section, in pti and pr 2 ■ We see a strong ppi dependence of the polarized 
contribution in combination with a large absolute size, which starts at 0% for pn = 1 GeV and goes up to 60% of 
the unpolarized for approaching 12 GeV, with p = 2m and Qq = 2 GeV. The results with p = mp have less 
polarization, with a maximal ratio reduced to about 10%. Some of the ppi dependence remains but most of it is at pxi 
values below 3 GeV. A precise measurement of this double differential cross section could be able to distinguish some 
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Figure 9. Normalized double differential cross section vs the transverse momentum of the two charm quarks in the extended 
region down to pn = 1 GeV, at p = 2m (left) and p = mr (right). The lower panels show the relative size of the polarized 
contribution. 


of the different scenarios, and either see first indications of or set first limits on the longitudinally polarized gluons in 
DPS. In particular if it is possible to extend the measured region down to lower transverse momenta. Investigating 
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Figure 10. Normalized cross section vs the azimuthal angle between the two charm quarks, at p = 2m, with kinematic region 
extended down to pn = 1 GeV. 
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the dependence of the relative size of the mixed unpolarized - linearly polarized gluons on the lower limit of transverse 
momenta we expect a rather large increase in the relative size (compared to the cross section contribution without 
azimuthal dependence). This is also visible in figure 10, where the relative size of the mixed contribution is increased by 
almost an order of magnitude when extending the kinematic region down to ppi = 1 GeV. The amplitude is still small 
compared to the angular modulation in the data, but it is another indication that allowing for a non-zero transverse 
momentum of the initial gluons, through for example NLO correction, could lead to significant enhancements. 
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A. Predictions at y/ s = 13 TeV 


In this section we show predictions at a hadronic center of mass energy of 13 TeV. The results are generally very 
similar to those at y/s = 7 TeV, and we will therefore keep the discussion rather brief. Figure 11 (first row) displays 
the normalized cross section as a function of the pt of one of the two quarks. The change in CM energy as compared 
to figure 3 flattens the cross section slightly and leads to a small decrease of the polarization. The second row of figure 
11 shows the cross section results as a function of the rapidity difference Ay. The change in CM energy has no visible 
impact on the shape of the cross section, and only leads to a small decrease of the polarized contribution. This small 
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Figure 11. Collisions at y/s = 13 TeV. Top: normalized cross section vs the transverse momentum of one of the charm quarks 
with y = 2m (left) and y = mr (right). The lower panels show the relative size R of the polarized contribution. Bottom: 
normalized cross section vs the rapidity difference Ay between the two charm quarks with y = 2 m (left) and y = mr (right) 
and relative size R of the polarized contribution. 
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decrease of the polarization as well as the small, if any, changes to the shape of the cross section is observed also for 
the dependence on M cc and A</>, as demonstrated by figure 12. 
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Figure 12. Collisions at s/s = 13 TeV. Top: normalized cross section vs the invariant mass of the two charm quarks with p = 2m 
(left) and /j, = tot (right). The lower panels show the relative size R of the polarized contribution. Bottom: normalized cross 
section vs the azimuthal angle between the two charm quarks A(j> between the two charm quarks with p = 2 m (left) and 
p = my (right) and relative size R of the polarized contribution. 


In the extended pxi region for the double differential cross section, the results have large polarization with a strong 
dependence on the transverse momentum at p = 2m and Q o = 2 GeV, as shown in figure 13. The contribution of the 
polarization decreases, as does the shape dependence, when going to /i = tut- 
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Figure 13. Collisions at y/s = 13 TeV. Normalized double differential cross section vs the pr of the two charm quarks, with 
p = 2 m (left) and p = mr (right), in the extended region down to pn = 1 GeV. 


VI. CONCLUSIONS 

We have investigated the effects of polarization in the double open-charm cross section, when the two charm quarks 
are produced in the kinematic region probed in the D°D° measurement by the LHCb Collaboration [12]. Polarization 
can give sizable effects on the magnitude of the cross section, reaching up above 50% of the unpolarized contribution 
in certain kinematic regions. The size strongly depends on the choices made when modeling the polarized double 
gluon distributions and on the large uncertainties for the single gluon distributions at the relevant low scales and 
small momentum fractions. We have presented the results obtained with y/s equal to both 7 and 13 TeV. The change 
of energy scale only has minor impact on the shape of the DPS cross section results as well as the relative size of the 
polarized contributions. 

The shape of the polarized contributions to the cross section are in most variables quite similar to the unpolarized 
results. In these cases it is difficult to disentangle the polarized contribution from other contributions in the DPS 
cross section, such as a single parton splitting and color interference contributions. We therefore identify variables 
and kinematic regions where the polarization does introduce some shape dependence. The most prominent shape 
dependence is found for the cross section double differential in the px of the two charm quarks, where the polarized 
contribution can vary with px from 0 up to 60% of the unpolarized. 

We compare the results of our calculation with the measurement of D°D° mesons by the LHCb [12]. For most 
distributions, the leading order calculation reproduces the experimental data rather well. The data cannot discriminate 
between the different models for the polarized DPDs as the polarization does not introduce any strong shape changes. 
The exception is the dependence on the azimuthal angle between the two mesons, which exhibits an approximate 
cos 2A (j> modulation. Polarized double parton scattering naturally produces such a modulation in the combination of 
linearly polarized and unpolarized gluons. However, the leading order DPS cross section for this term is too small 
to reproduce the modulation in the data. It is possible, however, that the size of this term changes drastically when 
including higher orders. Higher order effects for the process are expected to be large [17]. Such large NLO corrections, 
in combination with an expectation that the higher order corrections will lift the strong suppression of the mixed 
(unpolarized - linearly polarized) contribution present at tree level, can lead to a significant enhancement of the 
amplitude of the azimuthal modulation. Unfortunately, the theoretical formalism for the description of DPS needs to 
be further developed to reach a state where higher order effects can be systematically included. 

In the double differential cross section, looking at the px of both of the charm quarks, the longitudinal polarization 
can have a larger impact on both the size and shape. Measurements of this double differential cross section could 
therefore give first experimental indications of, or limits on, the effects of polarization in double parton scattering. 

We have used the homogeneous double DGLAP evolution equations, which do not include any single part splitting 
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term. The effect of the single parton splitting on the unpolarized DPS cross section was studied in [15]. Including it 
also for the polarized terms of the DPS cross section could further enhance the effects of polarization. In addition, we 
have employed an ansatz which splits the unpolarized gluon DPD in two single gluon PDFs and a factor depending 
only on the transverse distance between the two partons. This approach, common in DPS studies, is useful as a 
first approximation of the gluon DPD, but neglects several effects. These include correlations between kinematical 
variables and the color of the two gluons. The cross section ratios which we present, are likely to be more stable to 
such corrections than the absolute size of the cross section. However, further phenomenological as well as experimental 
studies are required to better constrain these effects. 
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